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§0 Abstract 

We establish a complete set of hydrodynamic equations for a spin-1 Bose-Einstein condensate 
(BEC), which are equivalent to the multi-component Gross-Pitaevskii equations and expressed in 



00 



terms of only observable physical quantities: the spin density and the nematic (or quadrupolar) 
tensor in addition to the density and the mass current that appear in the hydrodynamic description 
of a scalar BEC. The obtained hydrodynamic equations involve a generalized Mermin-Ho relation 
that is valid regardless of the spatiotemporal dependence of the spin polarization. Low-lying collec- 
tive modes for phonons and magnons are reproduced by linearizing the hydrodynamic equations. 
We also apply the single-mode approximation to the hydrodynamic equations and find a complete 
set of analytic solutions. 

PACS numbers: 03.75.Kk, 03.75.Mn, 05.30.Jp 
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I. INTRODUCTION 



The standard mean-field description of scalar Bose-Einstein condensates (BECs) is given 
by the Gross-Pitaevskii (GP) equation [l, 2] which has proved to be quite instrumental in 
describing various properties of the system such as collective modes (3 6 j and dynamical in- 



stabilities 



llj . The hydrodynamic theory provides an equivalent yet intuitively appealing 



description of the system, since the equations of motion are expressed in terms of physical 
quantities such as the density of particles and the mass current [12J, [l3|. The mean- field the- 
ory of spinor BECs has been developed based on multi-component GP equations by Ho [ljj] 
and Ohmi and Machida It is natural to ask what are the corresponding hydrodynamic 
equations and what physical quantities other than the density of particles and the mass 
current are needed to make a complete mean-field description of the spinor BECs. Spin do- 
main formation and texture in a quenched spin-1 BEC were observed 
high-resolution imaging technique for magnetization profiles for spinor BECs 
prompted the hydrodynamic description for ferromagnetic spinor BECs 



by in-situ and 



20 



m 
22 



19| . which 



However, 



the Berkeley experiment 



161 ] shows that the magnetization is not fully polarized over the 



entire condensate. Meanwhile, the Majorana representation has been employed to describe 
general spin states 23|, |24j. Barnett et. al. 23[ have developed the mean-field hydrodynamic 



equations that involve the Landau-Lifshitz equations for the spin-node vectors and reproduce 
collective excitations from the viewpoint of the point-group symmetry, while Lamacraft 24 1 
has derived the low-energy Lagrangian and obtained the spin-wave spectra. In this paper, we 
derive the most general mean-field hydrodynamic equations for spin-1 BECs that are equiv- 
alent to the multi-component GP equations and expressed in terms of observable quantities 
such as the magnetization vector and the nematic (or quadrupolar) tensor. 

We consider a BEC of spin-1 bosons which interact via an s-wave contact interaction. As 
shown Sec III II a general mean-field state of a spin-1 BEC is completely characterized by the 
spin and the nematic (or quadrupolar) tensor in addition to the density and mass current. 
The density p, the spin F, and the nematic tensor N^ u are defined as 



P = 



(1) 

(2) 
(3) 



where a field operator ift^ annihilates a boson with spin index /i, F is the vector of spin-1 
matrices, and is a rank-2 symmetric tensor defined by 

= \ (W XT (F v )rr, + • W 

Throughout this paper, we employ the Cartesian representation and assume that repeated 
indices are to be summed over the Cartesian coordinates (see Sec. [II]) . The nematic tensor, 
which was originally introduced as an order parameter of liquid crystals, has attracted 
renewed interest in the field of ultracold atomic e; probe of the nematic order 

and in experiments of spin-nematic squeezing in an SU(3) system {25], 26]. 
The second-quantized Hamiltonian of the system is given by 



H = T + U + H Z + V, 



(5) 



where T, U , Hz, and V represent the kinetic-energy operator, spin- independent one-body 
potential, Zeeman terms, and two-body interaction, which are given as follows: 

h 2 _ 



T = J drip 
U = 



2M 



Hz = J dr$l(-pn ■ F + q(n ■ F) 2 ) 
dr (c : p 2 : +c x : F 2 :V 



(6) 

(7) 
(8) 

V = J dr {c : f> z : +c, : F' z :J . (9) 

Here M is the mass of a boson, U(r) is a spin-independent potential such as an optical 
confinement trap, p and q are the coefficients of the linear and quadratic Zeeman energies, 
respectively, Co and c\ represent the spin-independent and spin-exchange interaction energies, 
and the normal order of operator A is denoted as : A :. 

The mean-field dynamics of the system described by Hamiltonian ([5]) is governed by the 
time-dependent multi-component GP equations: 



ih 



dt 



h 2 



2M 



V 2 + U(r) 



+ 



p{F M ) + q(N zz ) 



cqS^ + cif x {F x 



/IV 



(10) 



where ifj^ represents the order parameter of the condensate, is the Kronecker's delta, and 
p and fn represent the particle-number density and the spin density: 



P = 



(11) 



U = C(F,) vX Cx, (12) 
with being the normalized condensate wave function defined by 

C„ = MVP- (13) 

In the following sections, we will derive a set of hydrodynamic equations that are equivalent 
to Eq. (jTUI) . and demonstrate how to reproduce various properties of the spinor BEC from 
them. 

This paper is organized as follows. Section [TT] derives the hydrodynamic equations and a 
generalized Mermin-Ho relation of spin-1 BECs. Section ITTT1 shows the equivalence of the de- 
rived hydrodynamic equations to the multi-component GP equations. Section IIVI and [V] are 
devoted to applications of the obtained hydrodynamic equations: low-lying collective modes 
are reproduced in a physically transparent manner in Sec. IIVt and the single-mode approx- 
imation is employed to the hydrodynamic equations to obtain a complete set of analytic 
solutions in Sec. [V] Finally, Sec. |VT] summarizes the main results of this paper. 



II. HYDRODYNAMIC EQUATIONS OF SPIN-1 BECS 

We adopt the Cartesian basis =: C (/i = x,y,z) to express operators in terms of 

their matrix representations. Each element of C satisfies 

% l/i) = 0. (14) 
The spin matrices in the Cartesian representation are given by 

Wu\ = -i^A, (15) 

where e^ v \ is the completely antisymmetric unit tensor of rank three. Then, the matrix 
elements of the nematic tensor defined in Eq. (J4]) reduces to 

(Np V ) Xr) = S^Sxr, - - {5^u V + <^a<W • (16) 

The hydrodynamic equations for spin-1 BECs are written down in terms of density p, 
mass current v, spin density f^, and nematic tensor n^, where the mass current and the 
nematic tensor are defined as 

^2m[C(vc,)-(vc;)cJ, (17) 
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and 

^=a(^) Ar? c„ (is) 

respectively. These variables p, v, and n^ v are referred as the hydrodynamic variables in 
this paper. The spin and nematic currents are defined as: 

*V = ^W.a [C (VCa) - (VC) Ca] , (19) 

*V = ^Wa, [Ca (VC,) - (VCa) C„] • (20) 

We shall express them in terms of the hydrodynamic variables by making use of the following 
identity: 



Q&> = <V - »v + 2 e ^fx- ( 21 ) 



We derive the hydrodynamic equations by substituting the time-dependent multi- 
component GP equation into the time derivatives of the variables p, v, / M and n^. Firstly, 
the time derivative of the density p leads to the mass continuity equation, i.e., 

^ + V-pv = 0, (22) 

which takes the same form as the mass continuity equation for scalar BECs and ferromagnetic 
BECs. 

Secondly, we obtain the continuity equation for the spin density as follows: 

^A + V • pVf, = jfz^vp (pfv ~ 2qn zu ) . (23) 
The spin current can be expressed in terms of the hydrodynamic variables as follows: 



^/i/(V/x) +n vv {Vn Xv ) 



(24) 



where the first term on the right-hand side arises from the drift of atoms with the spin density 
/, the second and the last terms describe the spin currents driven by the spatial variations 
(i.e., textures) of spin and spin nematicity, respectively. The last term in Eq. (1241) vanishes 
in the case of a fully-polarized (or ferromagnetic) BEC. In the general spin-1 case, however, 
spinor properties are described not only by the spin density but also by the nematic tensor 
as detailed in the next section. We also note that Eq. (|23|) reduces to what was obtained in 



Rcf. 



3- 22 



22J in the limit of the fully-polarized state, i.e., \f(t)\ = 1. 



Thirdly, the continuity equation for the nematic tensor, which does not appear in the 
case of a ferromagnetic BEC, is given by 



+ V ■ pv^ = -p 



dt ' r h 



(25) 



with the nematic current v^ v given by 
h 



AM 



{e M Ar? [fx (V%) - (V/ A ) n uv \ + e vXrj [f x (Vn m ) - (V/ A ) n m ]} . (26) 



The last two terms on the right-hand side of Eq. (1251) act as spin torques on the nematic 
tensor, for they can be rewritten as 



Cl 

h 



P 2 ((/ x nu), + (/ x »„)„) , (27) 



where n M = {n^^n^y^n^) . The nematic current given in Eq. (126]) involves the texture of 
this torque force in addition to the mass current associated with spin nematicity. We will 
discuss the effect of the torque force in Sec. [V] 

Finally, the equation of motion for the mass current is given by: 



— h (ViS/i) Vj -Vi — 

dt J) 2M 2 1 Jp 



(28) 



+ [(Vin^) (Vjn^) - {ViVjUft,,)] 

= ~M~ K V * U W + Co + Cl/ " ( V ^)] • 

This equation may be regarded as the Euler equation for the spin-1 BEC. On the left- 
hand side, the first two terms represent the usual material derivative and the third one is 
the quantum-pressure term; the remaining terms show the contributions from the spin and 
nematic textures. The right-hand side shows the force terms arising from the gradients of 
the one-body potential, the particle density, and the spin density, respectively. 
The vorticity, which is the rotation of the mass current, is given by 

V x v = (V x + rv ( v x *V) • ( 29 ) 



Equation ( 129]) may be regarded as a generalized Mermin-Ho relation which is valid for an 
arbitrary state. This result shows that both spin and nematic currents contribute to the 
vorticity through their rotations. In the fully-polarized limit, Eq. ( 12 9 p reduces to: 



L (V x v. 



0. 



(30) 



Equation (130]) is another expression of the Mermin-Ho relation 271 ] : 

h 



V x v 



2M 



e^xU {WU x V/ A ) = 0, 



(31) 



which can be confirmed by substituting Eq. (|24j) and n^ u = (5^ u + f^f u )/2, which holds for 
fully-polarized BECs, into Eq. ([30]). 

For use in Sec. [V] we derive the expression of rewrite the energy functional in terms of 
hydrodynamic variables. We start with the mean- field Hamiltonian corresponding to Eq. §5§, 
and make the replacement Eq. ( fTBT) . Using Eq. ( 12T1) . we obtain 



Co 



Cl 



-Ppfz + qpn zz + -^p A + -^p 2 ff 



(32) 



Here, the terms in Eq. ( 132]) can also be derived from the symmetry argument similar to the 
Ginzburg-Landau free energy as follows. The system is rotationally invariant except for the 
linear and the quadratic Zeeman terms so that the energy functional should be written as 



E 



dx 



(33) 



where Xp, X/; Xn, Xint, and x' nt are functions of p. Here we assume that there are no terms 
higher than the second order in hydrodynamic variables and a term of the second order in 
the nematic tensor. Since 

1 



(34) 



a term proportional to r? is absorbed by Xint and x( n t- Thus, with proper identification 
of the coefficients, Eq. (1331) reduces to Eq. (13"2"]) . Since only the invariant properties in spin 
space are used in deriving Eq. (!33j) . it might be used for the finite-temperature theory of a 
spinor BEG. 



FIG. 1. Euler angles a, /3, and 7 specifying an arbitrary rotation of the condensate wave function 
in the real space. The Euler rotation of a spinor wave function at the coordinate Cfj,(R) is described 
by the matrix representation of the rotation operator U(a(R), (3(R), 'j(R)) given in Eq. (|36p . The 
vector / represents the spin vector. 

III. COMPLETENESS OF THE HYDRODYNAMIC EQUATIONS 

In the case of scalar BECs, the hydro-dynamic equations are written in terms of particle- 
number density p and mass current v, and describe the same dynamics as the GP equation. 
The scalar hydrodynamic equations are therefore complete in this sense. It is natural to ask 
whether the hydrodynamic equations for the spin-1 BECs obtained in the preceding section 
can describe the same complete mean-field dynamics as the multi-component GP equations. 
In this section, we answer this question in the affirmative. 

An arbitrary mean-field state of a spin-1 BEC can be obtained by an Euler rotation 
of a state which is, in general, partially polarized in the z direction. Thus, a condensate 
wave function is determined by six variables, which give the same degrees of freedom as 
the three-component condensate wave function. Those new variables involve the particle- 
number density p, the phase of the U(l) gauge 0, the Euler angles a, f3, and 7 defined in 
Fig. [H and the polarization parameter ^ of a ^-polarized normalized wave function 

^(a,/3, T ;0) = V%*(«>fr7;0) = ^[U{a, /?, y)]^ z (0). (35) 

Here, the matrix [W(ot, 7)] indicates the Cartesian representation of an Euler rotation 
operator given by 

U(a, (3, 7) = e-^e-^e"^ 7 , (36) 
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and (J[ z is written as 



= costf, & (i?) = zsintf, Cj*(0) = 0, 



(37) 



where the polarization parameter d parametrizes the three phases of the spin-1 BEC (n G Z): 







2n+l 



2 71 " 



7T 



(/^ = F 2 = 1, fully-polarized or ferromagnetic); 
(/^ = 0, non-polarized or polar); 



(38) 



otherwise (0 < f 2 < 1, partially polarized). 



We shall refer \U(a, (3, 7)]„„d' Z ('#) lT1 Eq- (03) as the spinor part of the condensate wave 
function. Then, all the hydrodynamic variables can be rewritten in terms of the six variables 
of the wave function. Here, we examine the spin density f^, the nematic tensor n^ u , and the 
mass current v, since the density p is just the squared magnitude of the wave function. 

First, we express the spin density and nematic tensor in terms of the Euler angles 
a, j3, and 7 and the polarization parameter d. The spin density is calculated from Eq. ( fl2l) 
as 

/ = sin2^ e/ , (39) 

where the unit vector ej points in the radial direction of the unit vectors of the spherical 
coordinates, i.e., 



^cos a sin (3^ 



sin a sin f3 

cos j3 



/ 'cos a cos @\ 



sin a cos 
— sin (3 



f ' — sincA 



V 



COS OL 





(40) 



In terms of these basis vectors, the nematic tensor given in Eq. ffl8|) is diagonalized as: 



n^u = X r e ril e ru (r is summed over r = 1, 2, 3), 



(41) 



with the eigenvalues and eigenvectors given by 



Xi = \ (1 - cos2tf) , A 2 = ^ (1 + cos2tf) , A 3 = 1, 



and 



(42) 



ei = cos ^yep + sin7e a , e 2 = — sinjep + cos 7e Q , e 3 = e^. (43) 

Here, the eigenvector e3 is parallel to the spin vector whose magnitude can be expressed 
in terms of the eigenvalues of the nematic tensor as 2\ / X 1 \ 2 . On the other hand, the 




FIG. 2. Cartesian representation of the condensate wave functions at position R plotted in terms of 
r{6, 4>) given in Eq. (|44p with its two axes e.\ and e.<i indicated by red arrows. These two axes are on 
the plane delimited by closed contour perpendicular to the unit eigenvector of the nematic tensor 
e3 at R. For the sake of simplicity, here we set the Euler angles a = and /3 = so that the unit 
eigenvector of the nematic tensor e% is parallel to the z-axis and r(8,4>) = |C^(0, 0, 7; i?)V^(0, <f>)\. 
The polarization parameter $ is taken to be 1? = 37r/16. Note that e.\ and e.2 coincide with the 
major and minor axes of the condensate wave function, respectively. 

remaining eigenvectors e\ and e 2 are parallel to two of the three principal axes of the 
Cartesian representation of the probability amplitude (see Fig. [2]) given by 

r(M) = IC>,/?,7;WM)l, (44) 

where V^s are the basis functions of the Cartesian representation which are defined in terms 
of the spherical harmonic functions Y"/™(#,(/>) as 

V x (6, 4) = -j= (-Yiie, 0) + Yr\e, 0)) , 

V y (9A) = -j= Wfafi + Yf 1 ^)) , (45) 

K(M) = n°(M). 

The axes of the spherical coordinate plot of r (6, (ft) can be written in terms of the eigenvectors 
of the nematic tensor e% and e 2 as 



e r = wA(i_A r )e P (r = l,2), (46) 

V o7T 

which are shown in Fig. [2j These relations between the eigenvectors and eigenvalues of 
the nematic tensor and the probability amplitudes are analogous to the classical electric 
quadrupole moment that is caused by the distribution of the charge. 

10 



The spin density and nematic tensor have nine independent components, despite the fact 
that the numbers of the variables in the spinor part of the wave function are four, i.e., a, 
(3, 7, and The degrees of freedom of the spin density and nematic tensor, however, are 
four which is the same as the spinor part of the wave function because of the following five 
constraints: 

= 2, (47) 

n„ufu = ffii (48) 
1 

The mass current v defined in Eq. fTl7|) can be expressed as 



det rv = -Jl (49) 



h 

v = — [( V0) - (Va) sin 2$ cos f3 - ( V7) sin 20] , (50) 

which can be rewritten in terms of the phase of U(l) gauge <j), the spin vector f^, and the 
eigenvectors of the nematic tensor e rfl (r = 1, 2, 3) as 

v = V(f)- 7^^v\ [f^rv (Ve rA )] . (51) 

Thus, we can determine up to a constant factor from v and vice versa, if we know / M and 
n^y. Here, we note that the rotation of Eq. f)50p gives another expression of the generalized 
Mermin-Ho relation. Other expressions related to the generalized Mermin-Ho relation have 



also been discussed in Refs. 



23 



28 



29). 



We have thus shown that the hydrodynamic variables have the same degrees of freedom 
as the multi-component condensed wave function for a spin-1 BEC, and that a set of the 
hydrodynamic equations in Eqs. fl22l . fT2"3l . fl25l) . fl28|) . and fl29l) are equivalent to the multi- 
component GP equations. 

IV. COLLECTIVE MODES 

In this section, low-lying collective modes are derived from the hydrodynamic equations 



in the absence of an external magnetic field in a manner similar to those for a scalar BEC 12 , 



13(] . We show that the collective modes obtained from the multi-component GP equations 



with the Bogoliubov approximation are fully reproduced in a physically transparent manner. 
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A. Linearization 



The hydro dynamic equations are linearized as follows. The particle-number density p, 
the spin density f^, and the nematic tensor n^ u are decomposed into their c- number parts 
and fluctuations from them as 

p{t, r) = p + Sp{t, r), f^(t, r) = f li + <J/ M (t, r), n^(t, r) = + Sn^t, r). (52) 

The mass current v is assumed to be the first order in the fluctuations. Then, the second 
derivatives of the fluctuations of p and / M with respect to time are linearized up to the first 
order in the fluctuations as follows: 

SHp * (v 2 V^) + 2^Sp(V» + ^/„(V^), (53) 



and 



dt 2 AM 2 v rj M v ' M 



^ = -4^ C™'J + & ^ ~ ® ^ (54) 



+ (2n„„ - S„„ - /„/„) (V 2 <S/„) 



where we use the identity 



iV M AU = N ia/ -± {5, v fl - fj„) , (55) 



and the constraints in Eqs. (147j) - (149j) . The linearized equation for the nematic tensor n^ v 
will be discussed in the next two subsections. 

Collective modes can be decomposed into the density mode and magnon modes, since 
in the linearized modes the first-order terms in fluctuations are taken into account. The 
magnon modes involve two types of modes, i.e., the deformation mode and rotational mode 
with respect to the spin density and the nematic tensor. The magnitude of the spin vector or 
the shape of the nematic tensor varies in the deformation mode, while the orientation of the 
spin vector and the unit eigenvectors of the nematic tensor rotate in the rotational mode. 
In the magnon-deformation mode, the polarization d changes to d + <50, which implies that 
the fluctuations Sf^ and SX r (r = 1, 2) are given from Eqs. fl39|) and fj4~2]) as: 

5f„ = 25$ cos 20, (56) 

5Ai, 2 = ±50 sin 20. (57) 
Those relations will be used in the following subsections. 

12 



B. Ferromagnetic phase 



The spin density accompanied by the nematic tensor rotates in the rotational mode that 
is excited from the ferromagnetic state. The spin density and nematic tensor satisfy the 
following conditions in this case: 

fi = hW li = 0, (58) 
I (<V + fM . (59) 



n 



ft i' 



2 



Therefore, the linearized equation of motion for the fluctuation in the spin density is given 
by 

^ = -iS* ■ (60 » 

Thus, the dispersion relation of the rotational mode is given by that of a free particle with 
mass M: 

hhj = ±e k , (61) 

where e k denotes the energy of a free particle with mass M, i.e., e k = h 2 k 2 /2M. 

The deformation mode in the ferromagnetic phase is the mode of the nematic tensor (or 
the spherical plot of the condensate wave function shown in Fig. since the fluctuation 
of the spin density is second order in the fluctuation of the nematic tensor according to 
Eqs. (1381) . (1391) and (HOI) . The equation of motion for the nematic tensor is linearized as 



dHn " U ~ ^ 'V 2 V 2 5rv) - (V 2 5tv) - (jf) 2 p 2 5n, v , (62) 



dt 2 AM 2 

which gives the dispersion relation of the deformation mode as 



hw = ± \sk — 2cip\ , (63) 



indicating that a finite gap 2c\p arises. This mode is interpreted as the mode of SF_ 



which is the longitudinal spin-wave mode 
m = -1 (30] . 



Q, 



15| . and describes the excitation from m — 1 to 



The density mode is obtained from the linearized equation for the density fluctuation as 

Hence, the dispersion relation of the density mode is given by 

hu = ±\fsk [e k + 2 (c + ci)p], (65) 
13 



which is a Goldstone mode and has the same expression as that for a scalar condensate except 
for the spin-exchange interaction energy c\. This shift is due to the spin-gauge symmetry of 
the ferromagnetic phase. 



C. Polar phase 

In the polar phase, the conditions for the spin density and nematic tensor obey 

U = 0. ( 66 ) 

and 

n^xnxr, = Ufa,. (67) 

Then, the linearized equation of motion for the fluctuation in the nematic tensor is obtained 
as: 

^ = -^(^H0 + «») 

Hence the dispersion relation is given by 

hw = ±\fe k (e k + 2cip), (69) 

which is also a Goldstone mode. 

The deformation mode can be understood as the vibrational mode of the spin density. 
In the polar phase, the spin density and nematic tensor satisfy 

U = 0, 5f ll = ±28tie fll , (70) 

and 

Ai, 2 = 1, or A 2 ,i = 0, 1, 5Ai = 5\ 2 = 0{5d 2 ), A 3 = 1. (71) 
Thus, the linearized equation is given by 

# = -4JF( v2v ^ + ¥^ v24 «- (72) 

Hence, the dispersion relation is obtained as 

hw = ±\/e k {e k + 2cip), (73) 
14 



which has the same dispersion relation as that of the rotational mode in the polar state. 
The magnetization oscillates at the frequency given in Eq. (]73p . which can be understood 
as a consequence of the spin-exchange interaction. 

The two degenerate modes in Eq. (1B"91 an d E q. (1731) have conventionally been interpreted 



as linearly-independent magnon modes 



14 



15 



30] . Our hydrodynamic description reveals 



that they actually describe the rotation of the nematic tensor and the amplitude oscillation 
of the spin, respectively. 

Finally, the linearized equation for the density mode is obtained by substituting / M = 
into Eq. (|53J: 

^ -3^>> + jH^> ■ < 74 > 

Hence, the dispersion relation is given by 



ftjjj = ±^e k (e k + 2c p). (75) 

The density mode from the polar phase is also a Goldstone mode and has the same form 
as that of a scalar BEC, because the mass current in the polar state is given by the same 
expression as that of a scalar BEC, i.e., v = (/l/M)V0. 



V. SINGLE-MODE APPROXIMATION 

As we have seen in Sec. [TTJ, the dynamics of a spin-1 BEC is described by a set of the 
hydrodynamic equations in Eqs. 022p . f[2"3"j) . f[2"5"j) . and (|28|) . In the sinigle-mode approxima- 
tion, we assume that the hydrodynamic variables are spatially uniform. The hydrodynamic 
equations can then be analytically solved. Firstly, the continuity equations for the mass, 
spin vector, and nematic tensor given in Eqs. fl22|) . fl23|) . and (1251) reduce to the following 
equations: 

p = const., (76) 

df 

-q^ = tz^c (pf u ~ 2qn zu ) (/i = x, y), (77) 
f z = m = const., (78) 



dt 



+ {tfj,Xrifxn U r) + ^vXrifx'n^) 



(79) 
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where c = Cip/h, p = p/ch, and q = q/ch. Secondly, according to Eq. ( )32|) . the energy 
functional with the spatially uniform variables is written as 



E< 



SMA 



Vcp 



Hn zz + - 2 {fl + f y ) 



+ const. 



(80) 



where V denotes the volume of the system and the term "const." represents the time 
independent part of the energy functional. Thus, we define the quantity 8 as 

1 



8 = qn zz + -{fl~ m 2 ) 



(81) 



which remains constant in time. Finally, we derive the differential equations for the zz- 
component of the nematic tensor n zz and the Euler angle a defined in Fig. [TJ The other 
components of the spin vector and nematic tensor can be calculated according to Eqs. (I76|) - 
f J79|) from these two variables. Here, the linear Zeeman energy p causes only the Larmor 
precession. Thus, we set p = hereafter, since the system with finite p in the rotating frame 
of the Larmor precession with the frequency ub = P is equivalent to the system with p = 0. 
Then, the differential equation for n zz is derived from Eqs. (176l) -( J79l) as 



dn Z ; 
dt 



Ac 2 

■2 



-m m. 



l) 2 + 2 (qn zz - 8) 



- ( -q + 1 J n z 



2 



(82) 



4c7(ra 22 ). 



Equation 
nuber m. 



is equivalent to the equation for the population of the magnetic quantum 



31 



state derived in Refs. 
calculating the time derivative of a = tan" 1 f y j f x as 

da 2qm {n zz — 1 



32(. The equation of motion for a is obtained by 



dt qn zz 
Substituting Eq. ( 182]) into Eq. fl83|) . we obtain 



8 



da 
dn 77 . 



rn 
c 



8/q - 1 



(83) 



(84) 



^ f{n zz ) {n zz - 8/q) y/ 'f(n zz ) 
Solving Eq. ( )84|) . we find that a is written in terms of the linear combination of the first-kind 
and third-kind elliptic functions. 

A complete set of hydrodynamic variables can be analytically derived from Eqs. ( 1821 ). 
(|84|) . and an initial condition. In this paper, we solve these hydrodynamic equations in two 
special cases for the sake of simplicity. One is the case q = 0, i.e., there is no magnetic field; 
the other is the case q ^ and m = 0, which describes the system with a transverse spin a 
magnetic field. 



1(3 



A. q = 



The spin vector /„ remains constant when q = 0, and S in Eq. flHTl) becomes 



The differential equation of n zz in Eq. ( 182]) reduces to 



(85) 



<9n 2 . 



4c 2 



n. 



/ 2 + m 2 \ 2 £ 2 (l-/ 2 ) 



2/ 2 



(86) 



and a = const. Equation (186]) has a unique solution for a given arbitrary initial condition. 
Here, we take the initial condition such that 



f(0) = A 





V / 



, W/u/(0) 



1 

|(1+ M ) 

o o U 



\ 



(87) 



wo)^ 



where the magnitude of the spin vector f = sin 2-^(0) and uq = cos2$(0) = y 1 — /q 
satisfies < $(0) < tt/4, and the initial Euler angles a(0), ,5(0), and 7(0) are set to be 
a(0) = 0, 0(0) = 7r/2, and 7(0) = 0. The solution of the hydrodynamic equations is given 
by 

\ 

I (1 + Mo cos 2ut) ^uosm2cut 
^0 |u sin2u;i \ (1 — u cos 2ut) ,j 



\ 



(88) 



with 



u = c\f\. 



(89) 



The obtained nematic tensor in Eq. (188]) is diagonalized as = Xre^e^ (the sum is taken 
over r = 1, 2, 3), where 



Ai = ^ (1 + w ) , A 2 = ^ (1 - u ) , A 3 = 1, 



(90) 



and 



I \ 

cos cut 
ysinwt / 



e 2 



' \ 

— sin ut 
\ cosut J 



e 3 





V / 



(91) 
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which implies that the quadrupole moment rotates about the spin vector at frequency u = 
c\f\, despite the fact that there is no external magnetic field. This rotation of the nematic 
tensor is caused by the source terms on the right-hand side of Eq. (I79p ; the spin vector acts 
like a torque force on the nematic tensor because of the spin-exchange interaction. 

This rotational mode of the nematic tensor is not the magnon collective excitation 
from the partially polarized state 30[, but the mode caused by the spin- mixing derived 



in Ref. 



331 ] . The magnon collective excitation from the partially polarized state asymptot- 



ically approaches the deformation mode excited from the ferromagnetic phase derived in 
Eq. (I63p if we take the limit q — > with c < and \q/p\ > 0, while keeping q/p constant. 

Finally, we note that this phenomenon is expected to be observed in the following ex- 
perimental situation. First, a spin-1 BEC is prepared in a partially polarized state, for 
instance, in a BEC that consists of 80% particles in the m — 1 state and 20% particles in 
the m = — 1 state. Then, we excite particles in the m — — 1 (1) to the m = 1 (—1) state by 
a two-photon Raman process involving absorption of er + (er_) light and emission of o~_ (<t + ) 
light. If the frequency difference between two laser beams matches the transition frequency 
between m = — 1 and m = 1, we can expect a resonance phenomenon at u = c\f\. 



B. q / 0, m = 

A physical solution always exists in the case q ^ and m — 0. All hydrodynamic variables 
can be derived from the initial condition and n zz (t), since the angle of the spin vector in 
the x-y plane, a, which obeys the differential equation given in Eq. (1831) . remains constant. 
Thus, the equation of motion relevant to the dynamics is the differential equation for n zz 
given in Eq. (182)) . which is simplified as 

= 8C ^ { Uzz ~ f ) ( Uzz ~ A+) ^ Uzz ~ A_) ' (92) 
where two of three zeros, i.e., A±, are given by 

± 2 

There always exists a physical solution of Eq. (192]) with respect to an arbitrary initial 
condition, which can be analytically shown in a straightforward manner. The zz-component 




(93) 
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of the nematic tensor n zz obtained from Eq. ( |92|) is expressed in terms of the first-kind 
elliptic functions with elliptic modulus k determined from q and the initial condition, which 
implies that the dynamics has the period 2T(k), where T(k) is the first-kind elliptic integral. 
Unlike the case of q — 0, the magnitudes of the spin vector |/| and the axes of the spherical 
plot of the condensate wave function |ei| and |ei| defined in Eq. (I46p (see also Fig. [2]) depend 
on time. The eigenvectors of the nematic tensor e\ and e.^ rotate due to the torque force 
caused by the spin-exchange interaction as calculated for the case of q = 0; however, they 
do not always undergo circular motion, that is to say, ^(t) does not necessarily range from 
—7i to 7T, as a consequence of the energy balance between the quadratic Zeeman energy q 
and the spin-exchange energy c. Equation (1921) can be solved analytically with the same 
initial condition as in the case of q = in Eq. f[57|) . The solution depends on the condition 
for q and uq and the complete solution is listed in Appendix A and some typical examples 
are illustrated in Fig. 3. 



VI. SUMMARY 

In this paper, we have derived a complete set of the hydrodynamic equations for a spin- 
1 BEC in an arbitrary state, which give a description equivalent to the multi-component 
GP equations. The hydrodynamic equations are self-contained within the hydrodynamic 
variables: the particle- number density, the mass current, the spin density, and the nematic 
tensor, as shown in Sec. [Ill The energy functional in terms of the hydrodynamic variables has 
also been derived. The obtained equations involve the continuity equation for the density, 
the spin density, and the nematic tensor, the equation of motion for the mass current, and a 
generalized Mermin-Ho relation for arbitrary spin-1 BECs, as discussed in Sec. IIII1 We have 
applied our hydrodynamic equations to some specific situations in Sees. [TV] and [V] In Sec. [TV] 
we have shown that our hydrodynamic equations reproduce the collective modes obtained 
by the multi-component GP equations [lj, [uj] . In Sec. |VJ the single- mode approximation 
have been applied to the continuity equations so that the spin density acts on the nematic 
tensor like a torque as a consequence of the spin-exchange interaction. Under a magnetic 



field, the torque force and the quadratic Zeeman effect balance and t 
can be described by the first-kind elliptic functions as derived in Ref. 
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re dynamics of BEC 
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32|. 
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Appendix A: Solutions of Eq. (I92|) 



(i) — 1 — u Q < q < —2u , < q < 1 — u : 

The spin vector is always parallel to the x-axis because of the initial condition in 
Eqs. (187)) and the constant a and written as 

fl\ 

f(t) = f dn(cf t,k) 





v / 



(Al) 



where, dn(x, k) is defined in terms of the first-kind elliptic function sn(x, k) as 

dn (x, k) = y/1 - k 2 sn 2 (x,t), (A2) 



and k is given by 



k = 



^q{2u Q + q) 
fo 



(A3) 



Next, the nematic tensor is written as 



\ 



(A4) 



1 

1 - n zz (t) n yz (t) 
\0 n yz (t) n zz (t) J 

since the nematic tensor, which is a symmetric tensor with the condition in Eq. (jUJ, 
is diagonalized as Eq. ( 14T)) with the eigenvectors in Eqs. (143)) with e 3 = (1, 0, 0) T . The 
components of the nematic tensor are obtained from Eq. (1§2")1 and the yz- component 
in Eq. (J72J) as 

(A5) 



n vz{t) = r ( 2u o + q) sn (cf t, k) cn(c/ i, k), 



n. 



,{t) = - [1 - u + (2u + q) sn 2 (cf t, k)] . 



(A6) 
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Finally, the Euler angle j(t), which represents the rotational angle of two axes of the 
nematic tensor e.\ and e.2 about the spin vector, is uniquely determined by requiring 
that it satisfy 



lit) = 2 COS 1 



and 



lit) 



-sin 



uq cn 2 — (uo + q) sn 2 
UpCn 2 + (u + g) 2 sn 2 

(2u + q) sncn 



(A7) 



(A8) 



'liocn 2 + (uq + g) 2 sn 2 

where "sn" and "cn" are the abbreviations of sn (cfot, k) and cn (cfot, k), respectively. 
Taking the limit q — > + , Eqs (I A 71) and ( 1A8I) reduce to 7(t) in Eq. ( 1891) . that is to say, 
lit) = cf t. 

q = 1 — u : 

The elliptic modulus k = 1 in the case of q = 1 — Uq, which implies that the first-elliptic 
function sn becomes the hyperbolic function tanh. The spin vector is given by 

fo 



^ ^ cosh cfot 

The nematic tensor is written as in Eq. (1A4I) with n yz (t) and n zz (t) given by 

1 tanhc/ot 

2 cosh cfot 



and 



«zz(t) = - [1 - «o + (1 + «o) tanh 2 c/ t] 



(A9) 



(A10) 



(All) 



respectively. A BEC asymptotically becomes to the polar state in the limit t — > oo: 



lim f{t) = 0, 

t— >oo 



^10 0^ 



lim Uuvit) 



(A12) 





\° l / 

which implies that the axis of the nematic tensor in the polar state is along the y-axis, 
i.e., the direction perpendicular to the magnetic field. 



q > 1 - u , q < -1 — u : 

The spin vector is calculated 



as 



/(*) = /ocn (cf t/k,k) , 



(A13) 
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where 



k 



/o 



(A14) 



(2u + q) 

The nematic tensor has the same form as in Eq. (IA4I) . and n yz (t) and n zz (t) are written 
as 



n yz (t) = -^-sn(cf t/k,k)dn(cf t/k,k), 



n zz (t) = ^ ( 1 - u Q + ^sn 2 (c f t/k, k) 



Jo 
Q 



(A15) 



(A16) 



The angle j(t), which is the rotation angle of the eigenvectors of the nematic tensor 
ei and e2, satisfies 



/ % 1 -i 
7W = 2 C0S 



u dn 2 — (u + q) k 2 sv? 



Updn 2 + (w + qY k 2 sn 2 



~\2, 



(A17) 



and 



/ N 1 • -1 

TW = 2 sm 



f 2 
-/o 



dn sn 



qk 



Mgdn 2 + (w + q) z k 2 sn 2 



~\2, 



(A18) 



where "dn" and "sn" represent the abbreviations of dn(cf t/k,k) and sn(cf t/k,k), 
respectively. 



-2u < q < 0: 

The spin vector is expressed as 



/(f) = /ond(c/ t/vT^P,*) 



where nd(x, k) = l/dn(x, k) and the elliptic modulus k is given by 

1 



k 



1 - -r 



f 2 



q{2u +q) 



In this case, we have 



(A19) 



(A20) 



(A21) 



with the equality holding if and only if q = —Uq. For the sake of simplicity, the argu- 



ments (cfot/yl — k 2 , k) in elliptic functions are not written explicitly in the following 
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part of this example. The nematic tensor is written as the expression given in Eq. ( 1A4|) 
with 



n yz (t) 



2\q\VT^W 



=cd sd, 



n zz {t) 



1 



f 2 k 2 

l-u + 4^sd 2 



(A22) 



(A23) 
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where cd = cn/dn and sd = sn/dn. The rotation angle ^(t) is determined from the 
following equations: 



, * 1 , 

lit) = 2 C0S 



Mocd^ - (u + q) (1 - A; 2 ) sd 2 
_y / n 2 cd 2 + (n + g) 2 (l-/e 2 )sd 2 



(A24) 



and 



7 (t) = -Bin - 



cd sd 



|g| VT^ ^2^2 + {uo + ~ )2 (1 _ k2) sd : 
except for the case of q = —uq, where ^(t) is given by the following equations: 

t(*) = ^ cos_1 M ' 

7 (0 = ^sin" 1 [/ sd] . 



(A25) 



(A26) 



The eigenvectors of the nematic tensor ei and e 2 trace closed curves in the case of 
— 2uo < q < —Mo as shown in (1) of Fig. El whereas they rotates about e 3 when 
— u < q < as an example shown in (3) of Fig. |3l 

q — -2u : 

This case is trivial, where the spin vector and the nematic tensor remain constant. 
q = —1 — u : 

A BEC asymptotically becomes the polar state in the limit t — > oo. The spin vector 
is obtained as 



f(t) 



,fo 



(A27) 



cosh cfot 

and the nematic tensor, which is expressed as Eq. ( 1A4I) . involves the component n yz {t) 
and n zz (t) given by 

tanh cf t 



n yz(t) = -- (1 - Mo) 

2 cosh cj t 



and 



n zz {t) 



1 - u 
2cosh 2 cf t ' 



(A28) 
(A29) 



23 



(a) q=-2/5<-u 



?=0 



iXK) r(; f .j+o.oi 




T=T(k<)l2 t=T(k<) MT(k<)/2 



r(t<) 27-(fc) 



(b) q=-U 



t=0 




o rtfe) 2r(t) 



7{£ c )/2 f=r(& c ) 

z z 
A A 



teT(k c )/2 

A~ 



(c) 9 -l/4>-w ^ "hT(k>)l2 t=T(k>) t^3T(k>)/2 




o r(*>) 2r(i>) 



FIG. 3. Time and c/ dependences of \f(t)\ (each leftmost panel) and snapshots of the condensate 
wave functions with their principal axes ex's indicated by red arrows for the case of (a) q = —2/5 < 
Uo, (b) q = —uq, and (c) q = —1/4 > — uq, where g's satisfy — 2uq < q < 0. Here, 19(0) is taken 
to be #(0) = 37r/16 and hence uo = cos37r/8 and fo = sin37r/8. The period in each case is given 
by 2T(k(q)), where T{k{q)) is the first-kind elliptic integral with the elliptic modulus k(q) given 
by Eq. (|A20l) . We abbreviate (a) fc< = fc(-2/5), (b) /c c = k(-u ), and (c) £;> = fc(-l/4), here. 
The track of ex, which is indicated by solid closed curves with arrows, undergoes a transition at 
a branch point (b): q = —uq. The principal axes swing around their initial position for the case 
of q < —uq as shown in (a), whereas they rotate about the spin vector (or the eigenvector of the 
nematic tensor 63) for the case of q > —u$ as shown in (c), because for q = —uq, the major and 
minor axes exchange each other at t = T{k c ), where the spin vector is fully polarized and the major 
and minor axes can take any mutually orthogonal directions. 
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respectively. In the limit t — > oo, the spin vector and nematic tensor grow asymptoti- 
cally to 



lim f(t) = 0, lim n^{t) 



t— >oo 



t— S-oo 



d o (A 

1 
\0 Oy 



(A30) 



where the axis of the nematic tensor in the polar state is along the z-axis, i.e., the 
direction of the magnetic field. 
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